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SUMMARY 


Anisotropic plate theory was used to calculate the anisotropic stress concentra- 
tion factors for a composite plate (AS/3501-5 graphite/epoxy composite, single ply 
or laminated) containing a circular hole. This composite material was used on the 
highly maneuverable aircraft technology (HiMAT) vehicle. It was found that the ani- 
sotropic stress concentration factor could be greater or less than 3 (the stress con- 
centration factor for isotropic materials), and that the locations of the maximum 
tangential stress points could shift with the change of fiber orientation with 
respect to the loading axis. The effect of hole size on the stress concentration 
factor was examined using the Point Stress Criterion and the Averaged Stress Cri- 
terion. The predicted stress concentration factors based on the two theories com- 
pared fairly well with the measured values for the hole size 0.3175 cm (1/8 in). It 
was also found that through the lamination process, the stress concentration factor 
could be reduced drastically, indicating an improvement in structural performance. 


INTRODUCTION 


When an infinite isotropic plate containing a circular hole is subjected to 
remote uniaxial tension, the tangential stress along the boundary of the circular 
hole will reach a value three times the remote tensile stress at two boundary points 
lying on the hole diameter perpendicular to the loading axis. Namely, the maximum 
tangential stress concentration factor is 3, which is independent of the hole size. 
For anisotropic materials, such as fiber-reinforced composite materials, the picture 
is entirely different. The value of the maximum tangential stress concentration fac- 
tor for a composite plate can be greater or less than 3, and the locations of the 
maximum stress points could shift depending on the loading direction and the fiber 
orientations . 

This report calculates the tangential stress distribution around a circular hole 
in a composite plate (single ply or laminated) and examines how the maximum stress 
points shift with the fiber orientations. In addtion, the dependency of the stress 
concentration factor on the hole size is examined. In the analysis, the composite 
system is treated as a continuous anisotropic plate having effective elastic proper- 
ties. The calculated stress concentration factors are then compared with experimen- 
tal data. 


NOMENCLATURE 


a Q small fixed distance ahead of the hole boundary used in the Average • 

Stress Criterion 

d Q characteristic distance ahead of the hole boundary used in the Point 

Stress Criterion 

E-) modulus of elasticity of single ply in axis-1 direction or modulus of 

elasticity of anisotropic plate in axis-1 direction 
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modulus of elasticity of single ply in axis-2 direction or modulus of 
elasticity of anisotropic plate in axis-2 direction 

modulus of elasticity of laminated composite plate in axis-1 direction 

modulus of elasticity of laminated composite plate in axis-2 direction 

modulus of elasticity of anisotropic plate in a direction 

modulus of elasticity of single ply in fiber direction 

modulus of elasticity of single ply transverse to fiber direction 

shear modulus of single ply associated with {1,2} system 

shear modulus of single ply associated with {L,T} system 

stress concentration factor 

stress concentration factor at a = ir/2 



coordinate axis in fiber direction 

number of composite layers 

radius of circular hole 

rectangular Cartesian coordinates 

coordinate axis transverse to fiber direction 

rectangular coordinate system 

angular cooordinate 

angle between axis 1 and axis L 

complex roots of the anisotropic plate characteristic equation 

Poisson's ratios of single ply with respect to {1,2} system 

Poisson's ratio of single ply with respect to {l,t} system 

Poisson's ratios of laminated composite plate with respect to {1,2} 
system 

R/(R + d Q ) 
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5 2 R/(R + a Q ) 

a a stress in a direction 

Of tensile strength of composite plate 

a remote tensile stress 

03 

o y stress in y direction 

COMPOSITE ELASTIC CONSTANTS 


Let axes { 1 , 2 } be the coordinate axes, and let axes (l,t) be the principal 
elastic or material axes of the single composite ply shown in figure 1. The ply- 
elastic constants {e-|, E 2 , G-| 2 / v 12> v 2ll with respect to the { 1 , 2 } system can be 
related to the material constants {e^, Eip, G LT , v^ip, ^TL 1 with respect to the {l,t} 
system through the following equations (refs* 1 and 2). 
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If the composite plate is made of N number of single plies with different fiber 
orientations, then by using the mixture rule, the engineering elastic constants {e-|, 
e 2» Gi 2 , Vi 2 , v 2ll f° r the composite plate can be written as 
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The HiMAT composite plate (A5/3501-5 graphite/epoxy composite) is made up of 34 plies 
(N = 34, ply thickness = 0.0133 cm (0.00525 in)) with the following fiber 
orientations: 

14 plies of 0 = +50° fiber orientation 

14 plies of 0 = -50° fiber orientation 

6 plies of 0 = +35° fiber orientation 

The ply engineering elastic constants with respect to the principal elastic axes 


{l, T } are given by 

E L = 137.90 GPa (20 x 10 6 lb/in 2 ) (11) 

E t = 10.27 GPa (1.49 x l<)6 lb/in 2 ) (12) 

G lt = 2.41 GPa (0.35 x 1 0 6 lb/in 2 ) (13) 

v LT = 0.3 (14) 

Using equations (11) to (14), the ply elastic constants with respect to axes {1,2} 
can be calculated from equations (1) to (5). For 0 = ±50° fiber orientations, 

Ei(±50°) = 7.45 GPa (1.0809 x 10^ lb/in 2 ) (15) 

E 2 (±50°) = 8.43 GPa (1.2237 x 10 6 lb/in 2 ) (16) 

G 12 (±50) = 8.46 GPa (1.2274 x 10® lb/in 2 ) (17) 

v 12 (±50) = 0.5676 (18) 

For 0 = +35° fiber orientation, 

Et^S") = 9.59 GPa (1.3904 x 10 6 lb/in 2 ) (19) 

E 2 (+35° ) * 7.40 GPa (1.0732 x 10 6 lb/in 2 ) (20) 

G 12 (+35°) = 6.91 GPa (1.0025 x 10 6 lb/in 2 ) (21) 

*Note that the mixture rule does not give the relationship v-j 2 e 2 = v 2 iEi» 
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V-j 2 ( +35° ) = 0.6671 


( 22 ) 


Using equations (15) to (18) and equations (19) to (22), the engineering elastic 
constants for the HiMAT composite plate with respect to {1,2} system, can be calcu- 
lated from equations (6) to (9). 


14 14 a. 

e i =-jj e i < +5o °> E i ( “ 50O) + ti E i (+3SO) 

= 7.83 GPa (1.1356 x 10 6 lb/in 2 ) 

E 2 = jf E 2 (+50°) + E 2 (-50°) + E 2 (+35°) 

= 8.25 GPa (1.1972 x 106 ib/in 2 ) 

Gl2 = G 12 (+50°) + ^ G 12 (-50°) + G 12 (+35°) 


(23) 


(24) 


= 8.19 GPa (1.1877 x 10& ib/in 2 ) 

-14 14 6 

v 12 = 34 v 12 (+50°) + — v 12 (-50°) + — v 12 (+35°) = 0.5851 


(25) 

(26) 


TANGENTIAL STRESS AROUND A CIRCULAR HOLE 


For an anisotropic plate containing a circular hole subjected to remote uniaxial 
tensile stress a , acting at an angle <(> with respect to the principal elastic axis 1 
of the plate (fig. 2), the tangential stress, a a (or tangential stress concentration 
factor, K = °a/ a ) along the circular hole boundary may be expressed as (ref. 3, 
p. 174) 


(J ^ H! ^ 

K = — = — ]l_-cos 2 <f> + (k+n) sin 2 <j)J k cos 2 a 
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- n (1 + k + n) sin (j) cos <j> sin a cos a| 
where E a is the modulus of elasticity in the a direction (fig. 2) given by 
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(30) 


n = -i(P-| + W 2 ) = yf(^ ~ v 12 ) + ^ 

where i = >^7, and hi and P 2 are the complex roots of the anisotropic plate charac- 
teristic equation 
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For isotropic materials k = 1 and n = 2, and the stress concentration factor K 
(eq. (27)) reduces to 


K 


3 a 

— =1-2 cos 2(a - <f>) 


(32) 


which gives K = -1 at ( a - <j>) = 0 or it, and K = 3 at (a - 4>) = - tt/2 


HOLE SIZE EFFECT 


Consider a circular hole of radius R in an infinite anisotropic plate as shown 
in figure 3. If the remote uniform stress a m is applied in the y-axis (or axis 1) 

direction, then the normal stress O y in the y-axis direction at a point on the x-axis 
(or axis 2) in front of the hole may be approximated by (ref. 4) 



At the hole boundary x = R , equation (33) gives the stress concentrations factor 
K f/2 


K ir/2 


Oy(R,0) 


oftt/2, 4>=0 


1 + n 


(34) 


which agrees with equation (27) by setting a = ir/2 and <j> = 0. Equation (34) gives a 
constant value of 1^/2 for the same material regardless of the hole size. As the 
hole size decreases, the stress concentration factor Ktj/ 2 decreases and finally 
approaches unity (that is, a plate without a hole). For composite materials, the 
hole-size effect on the stress concentration factor K-n/2 becomes significant when 
the hole diameter becomes less than 3.048 cm (1.2 in) (ref. 4). 

To account for the hole-size effect in equation (34), two failure criteria were 
advanced by Whitney and Nuismer (ref. 5) for composite materials. Two theories are 
described in the following sections. 
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Point Stress Criterion 


The point stress criterion assumes that the failure will occur when the stress 
0 y ( x, 0 ) at a certain small fixed distance (or characteristic distance) d Q ahead of 

the hole boundary first reaches the tensile strength Of of the material (or tensile 
strength of the plate without a hole) (Fig. 3). Namely, 


a y (x, 0 ) 


= Of 

x=R+d 0 


(35) 


( 1 ) 

Using this criterion and equation (33), the stress concentration factor, K 1T ^2/ 
can be written as (ref. 5) 


(1) 1 2 3 

Ktt/2 = 1 + 2 ?1 + 2 
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where 


Si 


R 

R + d r 


For a large hole (that is ?i+1), equation (36) gives 

= k tt/ 2 = 1 + n 


K ir/2 




As the hole size decreases (that is, £-|-*0), equation (36) reduces to 


K (1) 

k tt/2 


Sl-K) 


= 1 


which corresponds to the case of no hole. 


137) 


(38) 


(39) 


Average Stress Criterion 


The average stress criterion assumes that failure will occur when the average 
value of 0y(x,O) over some small fixed distance a 0 ahead of the hole boundary first 
reaches the tensile strength of the material (without a hole) (fig. 3). Namely, 



(x, 0 ) 


dx 


Of 


(40) 


Using equations (33) and (40), the stress concentration factor K-ti/2 can be written as 
(ref. 5) 
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where 
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R 


R + a r 


For a large hole (that is, £ 2 + 1 ) , equation (42) reduces to equation (34), 

Ktt/ 2 = 1 + n 


r (2) 

K ir/2 


c 2 ->-i 


And for the case without a hole (that is, £ 2 -H)), equation (42) gives 


k ( 2 ) 

K ir/2 


?2^0 


(41) 
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(43) 
Namely, 

(44) 


(45) 


The values of the characteristic distances d^, and a Q are determined by means of curve 
fitting of the experimental data obtained from tensile tests of rectangular specimens 
containing holes of different sizes. From literature surveys (refs. 4 and 5), the 
values of d,-, and a Q for composite materials are in the ranges of 

d Q » 0.076 to 0.127 cm (0.03 to 0.05 in) 
a Q ® 0.381 cm (0.15 in) 


EXPERIMENTS 


The effect of hole size on the stress concentration factors predicted from 
the two theories discussed previously were examined by performing simple coupon 
tests. The widths W of the two rectangular specimens were W = 3.81 cm (1.5 in) 
and W = 1.905 cm (0.75 in). Both of the specimens contain small circular holes of 
2R = 0.3175 cm (1/8 in). The specimens were tested in a tensile test machine up to 
failure under uniaxial tension in the axis-1 direction, and the remote stress o m at 

the time of failure was calculated for each specimen. One specimen without a hole 
was also loaded up to failure in a similar manner to obtain the tensile strength. Of, 

of the material. The stress concentration factor Ktt/ 2 for each specimen with a hole 
was calculated from 


a f 
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00 


Ktt/2 - 
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RESULTS 


By replacing {e-j , E 2 , G 12 , v 12^' respectively, with {e-j, E 2 , G 12 , v^ 2 } in 
equations (27) and (28), the tangential stresses a a around a circular hole in a 
laminated HiMAT composite plate were calculated for three loading cases: <)} = 0 

(loading in axis-1 direction), <j> = tt/ 4, and <j> = tt/ 2 (loading in axis-2 direction). 
Figure 4 shows the plots of o a for both the laminated HiMAT composite (eq. (27)) and 
isotropic materials (eq. (32)) when the plate is under uniaxial tension in the com- 
posite elastic axis-1 direction. The maximum stress concentration factor K for the 
laminated HiMAT composite reached the value 2.43 (less than 3) at f( ur locations 
(a = ± 65° and a = ± 115°) instead of two locations (a = - ir/2) for he isotropic 
case. When the loading axis is <j> = tt/ 4 oblique to the composite axis 1 (fig. 5), the 
stress concentration factor K reaches the value of 3.49 (greater than 3) at two loca- 
tions (a = 135° and a = -45°). When the loading axis is parallel to the composite 
elastic axis 2 (fig. 6), the stress concentration factor K reaches the peak value of 
2.44 (less than 3) at four locations (a = ±25° and a = ±155°). 

For comparison purposes, similar calculations were made for a single ply of the 
HiMAT composite using equation (27). When loading is along the fiber direction 
(axis 1, fig. 7), the stress concentration factor K reaches 10.20 at four locations 
(a = ±80° and a = ±100°). When the loading is 45° oblique to the fiber direction 
(fig. 8), K reaches as high as 27.78 at two points (a = 114° and a = -66°). When the 
loading direction is transverse to the fiber direction (fig. 9), the value of K 
reaches 6.81 at four points (a = ±55° and a = ±125°). Tables 1 and 2 summarize the 
above results for both laminated and single-ply HiMAT composites. Notice that the 
laminated composite has much lower stress concentration factors than the single ply, 
indicating the improvement of the structural performance through the lamination 
process . 


Figure 10 shows the plots of equations (36) and (42) for the laminated HiMAT com- 


posite n = / 2(zr~ - V 12 J + taking the characteristic distances d Q = 0.1016 cm 

\ \ e 2 / G 1 2 J 

(0.04 in), d Q = 0.127 cm (0.05 in), and a Q = 0.381 cm (0.15 in) (refs. 4 and 5). For 


a hole diameter of 0.3175 cm (1/8 in), the predicted values of K-rr/2 are summarized 
below. 


^O' 

cm (in) 


cm ( in ) 


0.0106 (0.04) 1.4367 

0.1270 (0.05) — 1.3152 

0.3810 (0.15) 1.3498 


The stress concentration factors Ktt/ 2 determined from the uniaxial tensile tests 
(loading in axis-1 direction) of the laminated HiMAT composite specimens are tabu- 
lated below. 
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Specimen 


W, 

cm (in) 


2R, 

cm ( in ) 


k tt/2 


3.810 (0.75) 
1.905 (1.5) 


0.3175 

0.3175 


( 1 / 8 ) 

( 1 / 8 ) 


1.333 

1.220 


The above two data points are plotted in figure 10 and compare fairly well with 
the predicted curves, especially with the curves d Q = 0.127 cm (0.05 in) and a Q 

0.381 cm (0.15 in ) . 


CONCULDING REMARKS 


Anisotropic plate theory was used to calculate the anisotropic stress con- 
centration factors for single-ply and laminated HiMAT (highly maneuverable aircraft 
technology) research aircraft composite plates, each of which contained a circular 
hole. The analysis showed that the anisotropic stress concentration factor could be 
greater or less than three for anisotropic materials, (three is the stress concentra- 
tion factor for isotropic materials) and the locations of the maximum tangential 
stress points could shift with the change of fiber orientation in relation to the 
loading axis. It was found that through the lamination process the stress concentra- 
tion factor could be reduced drastically, and thereby the structural performance 
could be improved. The effect of the hole size on the stress concentration factor 
was examined in the light of the Point Stress Criterion and the Averaged Stress Cri- 
terion. The stress concentration factors calculated from these theories agree fairly 
well with the measured values for the hole size 0.3175 cm (1/8 in). 


Ames Research Center 

Dryden Flight Research Facility 

National Aeronautics and Space Administration 

Edwards, California, February 21, 1984 
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TABLE 1. - STRESS CONCENTRATION FACTORS FOR A LAMINATED HiMAT COMPOSITE PLATE, 
34 PLIES ( [ ±50 ° ] i 4 , [-50°] 14 , [+35° ] g ) 


Orientation of loading axis, <f>.. 

0° 

45° 

0 

O 

<71 

Stress concentration factor, K 

2.43 

3.49 

2.44 

Locations of peak tensile stress, a.. 

±65°, ±115° 

+135°, -45° 

±25°, ±155° 


TABLE 2. - STRESS CONCENTRATION FACTORS FOR A SINGLE PLY OF HiMAT COMPOSITE 


Orientation of loading axis, $ 

0° 

45° 

O 

O 

OA 

Stress concentration factor, K 

10.20 

27.78 

6.81 

Locations of peak tensile stress, a.. 

±80°, ±100° 

+114°, -66° 

±55°, ±125° 




Figure 1. Single com- 
posite ply under combined 
loading. 


Figure 2. Tension at an angle to a 
principal elastic axis 1 of an ani- 
sotropic plate with a circular hole. 
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Figure 4. Tangential stress distribu~ 
tion around a circular hole in HiMAT com 
posite plate [ 14 * (+50° ) s , 14 x (-50° ) s 
6 x (+35° ) a ]. * = 0°. 














Figure 8. Tangential stress 
distribution around a circular 
hole in a single-ply composite 
plate. = 45° . 


Figure 7 . Tangential stress distri- 
bution around a circular hole in a 
single-ply composite plate. <| > = 0 ° . 
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